ROTATING WAVES IN THE THETA MODEL FOR A RING 
OF SYNAPTIC ALLY CONNECTED NEURONS 
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Abstract. We study rotating waves in the Theta model for a ring 
of synaptically-interacting neurons. We prove that when the neurons 
are oscillatory, at least one rotating wave always exists. In the case 
of excitable neurons, we prove that no travelling waves exist when the 
synaptic coupling is weak, and at least two rotating waves, a 'fast' one 
and a 'slow' one, exist when the synaptic coupling is sufficiently strong. 
We derive explicit upper and lower bounds for the 'critical' coupling 
strength as well as for the wave velocities. We also study the special 
case of uniform coupling, for which complete analytical results on the 
rotating waves can be achieved. 



1. INTRODUCTION 

In this work we study rotating waves in rings of neurons described by the 
Theta model. The Theta model OH |SJ H>] , which is derived as a canonical 
model for neurons near a 'saddle-node on a limit cycle' bifurcation, assumes 
the state of the neuron is given by an angle 9, with 9 = (21 + 1)tt, I E Z 
corresponding to the 'firing' state, and the dynamics described by 

(1) a j- t = 1 - cos(0) + (1 + cofl(0))G8 + /(<)), 

where I(t) represents the inputs to the neuron. When (5 < this represents an 
'excitable' neuron, which in the absence of external input (7 = 0) approaches a 
rest state, while if (3 > this represents an 'oscillatory' neuron which performs 
spontaneous oscillations in the absence of external input. 

A model of synaptically connected neurons on a continuous spacial domain 
f2 takes the form: 

(2) ^§^- = l-co S (e(x,t)) + (l + co S (d(x,t)))\l3+g [ J(x - y)s(y,t)dy 

at L Jn 

(3) + s(x, t) = P(9(x, i))(l - cs(x, t)), 
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where J is a positive function and P is defined by 

oo 

(4) P(0) = W - ( 2l + 

l — — oc 

Here s(x,t) (x G fi,t € M) measures the synaptic transmission from the neuron 
located at x, and according to 0,® it decays exponentially, except when 
the neuron fires (i.e. when 6(x,t) — (21 + 1)tt, I € Z), when it experiences 
a jump. J5J says that the neurons are modelled as Theta- neurons, where the 
input I(x,t) to the neuron at x, as in Q, is given by 

I(x,t)=g J(x - y)s(y,t)dy. 
Jn 

J(x — y) (here assumed to be positive) describes the relative strength of the 
synaptic coupling from the neuron at x to the neuron at y, while g > is a 
parameter measuring the overall coupling strength. 

The above model, in the case c > 0, is the one presented in In the 

case c = this model is the one presented in |H] (Remark 2) and [jj]. We 
always assume c > 0. 

When the geometry is linear, Q = R, it is natural to seek travelling waves 
of activity along the line in which each neuron makes one or more oscillations 
and then approaches rest. In |Sj it was proven that for sufficiently strong 
synaptic coupling g, at least two such waves, a slow and a fast one, exist, and 
also that they always involve each neuron firing more than one time before it 
approaches rest, while for sufficiently small g such waves do not exist. It was 
not determined how many times each neuron fires before coming to rest, and 
it may even be that each neuron fires infinitely many times. Some numerical 
results in the case of a one and a two-dimensional geometry were obtained in 


In this work we consider a different possibility for the spacial geometry: 
il = S , so the neurons are placed on a ring and our equations are © and 

(5) = h(9(x, t)) + gw(6(x, £)) j* J(x - y)s(y, t)dy. 
with 

(6) h{9) = 1 - cos(0) + /3(1 + cos(0)), w(9) = 1 + cos(0), 
where J is continuous, positive and periodic 

(7) J(x) > VieR, 

(8) J(x + 2tt) = J(x) \/x £ M, 
and the solutions satisfy the periodicity conditions 

(9) 9(x + 2tt, t) = 9(x, t) + 2vrm Vx,ieM 
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(10) s(x + 2n,t) = s(x,t) Vx,ieR 

The integer m (the 'winding number') is determined by the initial condition 
6(x,0), and will be preserved as long as the solution remains continuous. 

In this geometry, a different kind of wave of activity is possible: a wave 
that rotates around the ring repeatedly. Such waves, that is solutions of the 
form: 

(11) 8(x,t) = (/>(x + vt) 

(12) s(x,t) =r(x + vt) 

where v is the wave velocity, are the focus of our investigation. 

In section [5] we show that in the case that the winding number m = 0, 
there can exist only trivial rotating waves. Thus the interesting cases are 
when m > 0. Here we study the case m = 1, the case m > 1 being beyond 
our reach. Thus, this work concentrates on the first non-trivial case. 

Our central results about existence, nonexistence and multiplicity of ro- 
tating waves can be summarized as follows (see figures 1 1121 for the simplest 
diagrams consistent with these results): 

Theorem 1. Consider the equations J3J),|0) with conditions \Q), $10\) . and 
m = 1. 

(I) In the oscillatory case /3 > 0: for all g > there exists a rotating wave, 
with velocity going to +oo as g — > +oo. 

(II) In the excitable case (5 < 0: 

(i) For g > sufficiently small there exist no rotating waves. 

(ii) for g sufficiently large there exist at least two rotating waves, a 'fast' and a 
'slow' one, in the sense that their velocities approach +oo and 0, respectively, 
as g — ► +oo. 

Therefore our results bear resemblance to those obtained in [Hj for the case 
of a linear geometry. We note that although for the rotating waves found here 
each neuron fires infinitely many times, the reason for this is that it is re- 
excited each time, because of the periodic geometry. During each revolution 
of the rotating wave, each neuron fires once, so naively one could think that 
the analogous phenomenon in a linear geometry would be a travelling wave 
with each neuron firing once - but this was shown to be impossible in [5]. 
It is interesting to note that while in .Hj some restrictions were made on 
the coupling function J, like being decreasing with distance, here no such 
restrictions are imposed beyond Q , JHJ ■ We would expect however that some 
restriction would need to be imposed on J in order to obtain stability of 
travelling waves. The whole issue of stability remains quite open and awaits 
future investigation. In the case ft = M, both numerical evidence in [3 |Hj 
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and results obtained in other models indicate that the fast wave is stable 
while the slow wave is unstable, so we might conjecture that this is true for 
the case investigated here as well - at least under some natural assumptions 
on J. Some analytical progress on the stability question in the case £7 = K 
has recently been achieved in 

Let us note that the model considered here, in the case (3 < 0, describes 
waves in an excitable medium, about which an extensive literature exists (see 
|1U| and references therein). However, most models consider diffusive rather 
than synaptic coupling. In the case of the Theta model on a ring, with diffusive 
coupling, and m = 1, it is proven in 0] that a rotating wave exists regardless 
of the strength of coupling (i.e. the diffusion coefficient), so that our results 
highlight the difference between diffusive and synaptic coupling. 

In section we reduce the study of rotating waves to the investigation 
of the zeroes of a function of one variable. In section 0] we investigate the 
special case in which the coupling is uniform (J(x) is a constant function), 
which, although artificial from a biological point of view, allows us to obtain 
closed analytic expressions for the wave-velocity vs. coupling-strength curves 
in an elementary fashion. We can thus gain some intuition for the general 
case, and obtain information which is unavailable in the case of general J, 
like precise multiplicity results. It is interesting to investigate to what extent 
the more precise results obtained in the uniform-coupling case extend to the 
general case, and we shall indicate several questions, which remain open, in 
this direction. In section we turn to the case of general coupling functions 
J, and prove the results of theorem ^ above, obtaining also some quantitative 
estimates: lower and upper bounds for the critical values of synaptic coupling 
coupling strength <?, as well as for the wave velocities. 

2. PRELIMINARIES 

We begin with an elementary calculus lemma which is useful in several of 
our arguments below. 

Lemma 2. Let f : ffi — > R be a differentiable function, and let 6, c £ R, b =/= 0. 
be constants such that we have the following property: 

(13) f(z) = c =>f( z )=b. 

Then the equation f(z) = c has at most one solution. 

PROOF: Assume by way of contradiction that the equation f(z) = c has at 
least two solutions zq < z\ . Define S C M by 

S={z>z \f(z) = c}. 

S is nonempty because z\ € S. Let z = inf S. By continuity of / we have 
f(s) — c - We have either z>z orz_=z , and we shall show that both of 
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these possibilities lead to contradictions. If z > Zq, then by l|13|l we have 

/(*>) = /GO = c 
sign(f'(z Q )) = sign(f'(z)) = sign(b) 

so we conclude that there exists Z2 G ( z o, z) with f(zi) = c, contradicting the 
definition of z. If z — zq then zo is a limit-point of 5, which implies that 
f'(zo) = 0, contradicting (|13|) . These contradictions conclude our proof. 

Turning now to our investigation, we note a few properties of the functions 
h(9) and w(9) defined by Jfjjl which will be used often in our arguments: 

(14) h((2l + 1)tt) = 2 VZeZ, 

(15) /i(2Ztt) = 2/3, VZ G Z, 

(16) w({2l + 1)tt) =0 VleZ, 

(17) w(22tt) =2 VleZ, 

Plugging (|11(I . H12|I into ©, and setting z = x + vt we obtain the 
following equations for (j)(z), r(z): 

(18) v<f>'(z) = h(<t>(z)) + gw{<t>{z)) f J(z ~ y)r(y)dy, 



(19) vr'(z) + r(z) = P{<f>{z))(l - cr(z)). 

In order to satisfy the boundary conditions l|9")l. (|10|l . <f> and r have to satisfy 

(20) (j)(z + 2tt) = cj)(z) + 2irm Vz G R, 

(21) r(z + 2tt) = r(z) Vz G M. 

Let us first dispose of the case of zero-velocity waves, v ~ 0. We get the 
equations 

(22) h(^z))+gw(^z)) f J(z-y)r(y)dy = 0, 

J — 7T 

(23) r{z) = P{(j){z)){l- cr{z)). 

If there exists some zo G K with </>(zo) = (21 + l)ir, I G Z, then, substituting 
z — zq into ill'L'l) and using (|14fl . (|16|) . we obtain 2 = 0, a contradiction. Hence 
we must have 

(24) <j)(z) ± (21 + 1)tt Vz G E, I G Z 
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which implies that P(<fi(z)) = 0, so that (|23|l gives r(z) = 0, and i12_?l) reduces 
to h(<f>(z)) = 0, and thus </>(z) is a constant function, the constant being a 
root of h(9). This implies, first of all, that the winding number m is 0, since 
a constant cj>(z) cannot satisfy (|20() otherwise. In addition the function h(9) 
must vanish somewhere, which is equivalent to the condition /3 < 0. We have 
thus proven 

Lemma 3. Zero-velocity waves exist if and only if m = and (3 < 0, and in 
this case they are just the stationary solutions 

r(z) = 

4>(z) = icos" 1 (j^y) +27rfc, fceZ. 

We will now show that the trivial "waves" of lemma |3| are the only ones 
that occur for m — 0. 

Lemma 4. Assume m = 0. 

(i) If (3 > i/iere are no rotating waves. 

(ii) If [3 < i/ie onZj/ rotating waves are those given by lemma\^ 

PROOF: Assume (v,(j>(z),r(z)) is a solution of JTSJ|,|Q5J| satisfying with 
m = 0, i.e. 

(25) <j){z + 2tt) = <p(z) Vz € M, 

and J2J. We also assume w ^ 0, otherwise we are back to lemma [3] We shall 
prove below that <j)(z) must satisfy (|24|l . and hence that P(4>(z)) = 0, so that 
by {IHl),lGU we have r(z) = 0, so that {THJl reduces to 

(26) «(/)'(z) = h(4>(z)). 

Since w ^ 0, if /i(6') has no roots ((3 > 0), l|26ll has no solutions satisfying (|25l) . 
If /i(6*) does have roots (/3 < 0) then the only solutions of (|2^|l satisfying (|2"5Tl 
are constant functions, the constant being a root of h(8), and we are back 
to the same solutions given in lemma [21 which indeed can be considered as 
rotating waves with arbitrary velocity. 

It remains then to prove that (|24|l must hold. Assume by way of contra- 
diction that <f>(zo) = (21 + 1)tt for some integer I. By JT1J , (JTBJ) , HHJ , and the 
assumption v =/= 0, we have 

cp(z) = (21 + 1)tt (f)'(z) = -. 

v 

Thus the assumptions of lemma with f — ip, c — (21 + 1)tt, b = | , are 
satisfied, and we conclude that the equation <fi(z) = (21 + l)ir has at most one 
solution, contradicting the fact that, by l|25|l . we have <f>(zo) = <fi(zo + 2ir) = 
(2Z + 1)tt. 
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Having found all possible rotating waves in the case m = 0, we can now 
turn to the case m > 0. In fact, as was mentioned in the introduction, we 
shall treat the case m = 1, the cases m > 1 being harder. By lemma [3] we 
know that there are no zero- velocity waves, so we can assume v ^ and define 

(27) A = i, 

v 

so that our equations for the rotating waves can be rewritten 

(28) <f>'{z) = \h{ct>{z)) + \gw{cj>{z)) f J(z - y)r(y)dy, 

J —TV 

(29) r'{z) + Xr(z) = XP(cj)(z))(l - cr(z)), 
with periodic conditions 

(30) 4>{z + 2rr) = cf)(z) + 2rr Mz e K, 

(31) r(z + 2tt) = r(z) Vz E M. 

3. Reduction to a one-dimensional equation 

We study the equations l|28|l . (|29|l for (A, 4>(z),r(z)) with periodic conditions 
(I30|l . (|31|) . We will derive a scalar equation (see (|52|l below) so that rotating 
waves are in one-to-one correspondence with solutions of that equation. 

We note first that, since by l|31(l we have </>(K) = K, and since any rotating 
wave generates a family of other rotating waves by translations, we may, 
without loss of generality, fix 

(32) 0(0) = 7T. 

The following lemma shows that for a rotating wave (in the case m = 1) 
there is at any specific time a unique neuron on the ring which is firing. This 
fact is very important for our analysis. 

Lemma 5. Assume (\,(j>,r) satisfy with conditions if%7J)) . ifff?)) . p%)) . 

Then 

(33) z e (0, 2tt) tt < <p(z) < 3tt 



(34) Z E (-27T, 0) -TT < (p(z) < 77 

PROOF: (£21 follows from (J^SJl by jgDjl. To prove (j22l, we first note that 
certainly A ^ 0, since A = and imply ^(z) is constant, contradicting 

We note the key fact that, by and l(TB|) . 

(35) 0(z) = (2^ + 1)tt, / e Z ^ 0'(z) = 2A. 
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By lemma (I35JI implies that the equation <f>(z ) = 21 + 1 has at most one 
solution for each I G Z. In particular, since 0(0) = tt, </>(27r) = 37r we have 
4>{z) ^ 7r, 37r for z G (0, 27r), and by continuity of <fi{z) this implies l|33|l . 

Let us note that if we knew that for rotating waves the function 4>(z) must 
be monotone, then lemma would follow immediately from 13211 . 

Question 6. Is it true in general that rotating wave solutions are monotone 
(for m^O)? 

Lemma 7. Assume (A, (f>,r) satisfy with conditions HI)) . Hi)) . Uity) . 

Then A > 0. 

In other words, v > for all rotating waves with m = 1, so the waves 
rotate clockwise. Of course in the symmetric case m = — 1 the waves will 
rotate counter-clockwise. 

PROOF: By (J2U) and 051) we have 0'(O) = 2A. We have already noted that 
A 7^ 0. If A were negative, then (j) would be decreasing near z = 0, so for small 
z > we would have (f>(z) < 0, contradicting (|33J) . 

Our next step is to solve (f2"§|l . (f3TJl for r(z), in terms of 4>(z). We will use 
the following important consequence of lemma 

Lemma 8. 

P{4>{z))\ { -2^) = ^S(Z) 

PROOF: By lemma|5]we have 

P{<f>{z))\ { ^, 2v) =8{cl>{z)-'K)), 

so we will show that 

(36) 5{<t>(z) - tt)) = ±.6(z). 

Let x( z ) G Co°(10 be a test function. Using lemma again we have 

(37) / x(u)(5(</)(it) — n)du = / x(u)5{4>{u) — Tr)du, 

J —TV J— € 

where e > is arbitrary. In particular, since 4>'(0) = 2A > 0, we may choose 
e > sufficiently small so that (f>'(z) > for z E (— e, e), so that we can make 
a change of variables ip — <fi(u), obtaining 

_ fan 
X(u)5(<j)(u) - ir)du = j x(0 _1 (^))^(^ ~ ^)-777— i 



(38) 



x(o) x(o) _ x(o) 



<t>>(<f>- x {-K)) <£'(0) 2A 
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This proves (|36|l . completing the proof of the lemma. 

By lemmalHlwe can rewrite equation l|29[) on the interval (— 2n, 2tt) as 

(39) r'(z) + (\ + ^S(z))r(z) = -\s(z), 
The solution of which is given by 

(40) r(z) = (^H{z) + r(-7r)e-" A ) e -(A*+*ff«) , 

where H is the Heaviside function: H(z) — for z < 0, H(z) — 1 for z > 0. 
Substituting z = ir into (|40|l and using l|31|) . we obtain an equation for r(— 7r) 
whose solution is 

r(-7r) = i(e' rA+ *- e - w V. 

and substituting this back into H4()|l . we obtain that the solution of H29(I . H31|I 
which we denote by r\{z) in order to emphasize the dependence on the pa- 
rameter A, is given on the interval (— 27r, 2tt) by 

(41) r A (z) = i e - (Az+ 5 ff ^))[i7(z) + (e 2 " A+ t-l)- 1 ] < \z\ < 2ir. 

We note that, for general z E R, r\(z) is given as the 27r-periodic extension 
of the function defined by (|41|l from [— 7r, it] to the whole real line. 

The following result, which can be computed from l|41|l . will be needed 
later 

(42) J r x (u)du = ^Pc(A), 
where 

2ttA _ -i 

Pc(A) 



e 27rA+§ _ 

We note that 

(43) p (A) = 1, 

A fact that considerably simplifies the formulas in the case c = 0. 
The rotating waves correspond to solutions (A, </>) of the equation 

(44) <p'(z) = \h(<f>(z))+\gw(ct>(z)) f J(z - y)r x (y)dy, 

with 4>(z) satisfying l|32|) and 

(45) (/.(tt) = 0(-tt) + 2tt. 
To simplify notation, we define 

(46) R x {z) = f J{z-y)r x {y)dy, 
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so that jilt is rewritten as 

(47) <t/(z) = \h((f>(z)) + XgRx(z)w((f>(z)). 

We note that l|47|l is a nonautonomous differential equation for <fi(z), and 
since the nonlinearities are bounded and Lipschitzian, the initial value problem 
(I47|l . (|32() has a unique solution, which we denote by <j>\. 

Rotating waves thus correspond to solutions A > of the equation 

(48) 4> x (tt) - 4> x {-tz) = 2tt. 
Rewriting (|47|l and 150(1 we have 

(49) <p' x (z) = + XgR x (z)w(4> x {z)), 

(50) A (O) = 7T, 

and defining 

(51) *(A) = -^a(tt) - 0a(-7t)), 

we obtain that rotating waves correspond to solutions A > of the equation 

(52) #(A) = 1. 

4. The case of uniform coupling 

Assuming that the coupling is J = 1 we shall be able to solve for the 
rotating waves explicitly. In this case we have, from 1(46(1 . 1(42(1 

r i 

R\{z) = J r x {y)dy = —p c (\), 
so that reduces to 

(53) <j>' x {z) = \h{4>x{z)) + 9 -p c {\)w{4> x {z)). 

The fact that l(53|) is an autonomous equation is what makes the treatment of 
the case J constant much simpler. Indeed, assume that 1(52(1 holds, so that 

(54) <£ A (7r)-0 A (-7r) = 27r. 

Then we have, using 1(53(1 . making a change of variables ip — A (z), an d using 
_ J_ /" r <t>'{z)dz 



(55) = 



2tt ./_„ \h(4>{z)) + £p c (\)w(<t>(z)) 

i /-^m dip i r 



2^ J M -ir) Xh(<p) + |p c (A)u)(v) 2^ AfcfoO + f Pc(A)w(^) 
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Substituting the explicit expressions for h and w from ©, and using the 
formula 

(155)) becomes 

(57) 1 ' 



so that rotating waves correspond to solutions of (|57() . with their velocities 
given by v = j. We can rewrite (IBTjl as 

(58) fcA*)=9, A>0 

where 

< 59 > M*> = 

In the following lemma we collect some properties of the functions / C ,/3(A), 
which are obtained by elementary calculus: 

Lemma 9. (i) When (3 < 0, f c g is positive and convex on (0,oo), and 
(60) jim / CiJ a(A) = oo, 

A — >Q 



(61) lim / C)1 a(A) = oo. 

A — >00 

(ii) When >0, f c p is decreasing on (0, oo), and ]60\) holds. If (3 > it has 
a zero at X — ^rp, if ft — it is positive on (0, oo) and linix— »oo / C) o(A) = 0. 

From lemma El we conclude that when (3 < l|58|l has exactly two solutions 
if g > 0(c, /?), where 

(62) fi(c,/3) =min/ Ci(3 (A), 

A^>0 

which we will denote by 

Ac^Cfif) < A ClJ a(ff), 

no solution if 5 < f2(c, /3), and a unique solution when g = f2(c, /?). 

When j3 > 0, part (ii) of lemma implies that (|58() has a unique solution 
for any g > 0, which we denote by A C)l a((7). 

An elementary asymptotic analysis of the equation (|58Jl yields 

Lemma 10. (i) When (3 < we have the following asymptotics as g — ► 00 

g ^ / 1 \ 

(63) A C:/3 (sO = ■^jj-S' + ^-J as.g-^00. 
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For X c pig), in case c > we have 

(64) X c ^(g) = y^(e £ 2 -l)^ + o{^) as g -> oo, 
while in case c = we have 

(65) \p{g) = \- + °{^) as g^oo. 

(ii) For /? > 0, ifte asymptotics of X c ,p(g) as g — ► oo are i/ie same as i/iose of 
Ac 3(3); ffiuen m |6'^| ) /or c > and |6'5|) for c = 0. 

We thus obtain 

Theorem 11. When J = 1; 
(TJ 7n i/ie excitable case (3 < 0: 

(%) 1/5 > f2(g,c) £/iere ea;zs£ two rotating waves with velocities given by 

(66) w c ,/3(.a) = = 1 , r , «c/3(.g) = . \ r , 
and we /iaue, /or i/ie sZow wave 

c 1 



(67) v c ^(g) = 2\P\ei- + 0[-) as g -> 00, 

/or i/ie /asi wave w/ien c > 0: 



(68) ^(ff) =2J-^—. -y/g + 0(l) asg^ao. 

V ea — 1 

wMe /or t/ie /asi wave w/ien c = 

(69) «o,/j(flO = 2 -? + Of-") as g^ 00. 

^g J 

(ii) If g = f2(c, /?) i/iere exists a unique rotating wave with velocity 

(70) v = v cS {g)=v c ^(g). 
(Hi) If g < f2(c, /3) i/iere exist no rotating waves. 

(II) When (3 > 0, there exists a unique rotating wave for any g > 0, whose 
velocity is given by 

(71) VcAd) = x 1 / v 

and /or Zarge g ii /ias i/ie same asymptotics as in jo'&)) . fo'^j) m i/ie cases c > 0. 
c = 0, respectively. 

In the excitable case we thus have two rotating waves born at a supercritical 
saddle-node bifurcation as the coupling strength g crosses Sl(c,{3). 
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Figure 1. Velocity of waves (v) vs. coupling strength (g) 
for the case J = 1, c = 0, (3 = —0.5. 



We now note that in the special case c = (the model introduced in 
we can obtain more explicit expressions. Using (|43|l we have 

/MA) = 2A 

The minimum in (|62|) can now be computed explicitly, and we obtain, when 
P <0, 



n(0,j9)=2 

We can also solve l|58|l explicitly, and obtain the velocities of the rotating 
waves. When (3 < 0, g > ft (0,(3) 



VjoAq) =9- vV + 4/3, v Ot0 (g) =g+ ^g 2 + 4/3. 
When (3 > 0, for all g > 



fo.Ms) = vV+4/3 + .9- 



Figures 1 1121 show the wave-velocity vs. coupling strength diagrams for 
the rotating waves when J = 1, c = 0, in an excitable ((3 = —0.5) and an 
oscillatory ((3 — 0.5) case. In figures Efl4l we change c to c = 1. 
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FIGURE 2. Velocity of waves (v) vs. coupling strength (g) 
for the case J = 1, c = 0, (3 = 0.5. 
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FIGURE 3. Velocity of waves (v) vs. coupling strength (g) 
for the case J = 1, c = 1, j3 = —0.5. 
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Figure 4. Velocity of waves (u) vs. coupling strength (g) 
for the case J = 1, c = 1, j3 = 0.5. 



5. The general case 

We now return to the case when J is a general continuous positive 2ir- 
periodic function, and prove that several of the results about rotating waves 
obtained above for the special case J = 1 remain valid, though the proofs are 
necessarily less direct. 



Lemma 12. 



lira *(A) = 0. 
x^o 



PROOF: We shall prove that 

(72) <j>x{z) =tt + 0(A) os A^O 

uniformly in z G [— tt,tt]. The lemma follows immediately from this and from 
When c > 0, the claim Q72|) is immediate, since, using (14111 . 

lim rx(z) = i e -^ H(z) [(e2 - l)- 1 + H(z)] , < \z\ < 2tt, 

so that 

(73) XRx(z) = 0(A) as A -> 0, 
which implies that the right-hand side of l|4T|l is O(X). 
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For c = 0, r\(z) becomes singular as A — > 0, so we need a more refined 
argument. For Xr\(z) we have 

(74) Xr x (z) = — + 0(X) as X -» 0, 

47T 

uniformly in z € [— 7T, 7r], hence 

(75) Afl A (*0 = ij + 0(A) as A^O, 
where 



- i r 

J = — / J(x)dx. 



The standard theorems on dependence of solutions of initial- value problems 
on parameters hence imply 

<j)\(z) = cj) (z) + 0(A) as X -> 0, 

uniformly in z G [— tt, it], where 4>o satisfies 

(76) <&{z) = ^w(Mz)) 

and ^o(0) = 7r. Since, by Ijltil) . the constant function n is a solution to 
this initial-value problem, the uniqueness theorem for initial-value problems 
implies that 4>o(z) = n. 

The following bounds, which follow immediately from l|46|l and (|42|l . will 
be useful: 

Lemma 13. (i) We have, for all z E [—ir,ir]: 

^-p c (X)minJ{x) < R x {z) < ^rPc{X) max J(x). 

ZA x£R ZA x£K 

(ii) If J is not a constant function, the inequalities of (i) are strict. 
Lemma 14. In the excitable case (3 < 0, we have 

— ^— > -^-max J(x) => *(A) < 1, 

PROOF: In the case of constant J the result can be proven by direct computa- 
tion, so we now assume J is not a constant function. We will show that when 

(77) — > -^-maxJ(i) 

V ; Pc{X) ~ 2|/3| xeR 1 ' 

we have 

(78) <j)\(z) < 2tt Vz e [0,7r], 
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(79) (t>x(z) >0 Vze [-77,0]. 

Together with 1)5 l|l. these imply the result of our lemma. To prove our claim 
we note that, using l)49)) . l)15)) . ()17)) . P ar t ( u ) °f lemma IT51 (which is why we 
need the assumption that J is non-constant) and (|77(l 

(f)x(z) = or 2tt => <f>' x (z) = 2X/3 + 2\gR x (z) 

(80) < 2A/3 + 5 p c (A)max J(a;) < 2A/3 + 2A|/3| = 0. 

We now show that (HQ implies (J7HJ). If JZSJ fails to hold, then we set 

zq = min {z G [0,7r] I </>a(-z) = 27r}. 

This number is well-defined by continuity and by the fact that <p\ (0) = 77, 
which implies also that zq > 0. By (|80)l we have 4>'\{zq) < 0, but this implies 
that (j)\( z ) is decreasing in a neighborhood of z$ , and in particular that there 
exist z € (0, zq) satisfying <j)\ (z) = 277. But this contradicts the definition of 
zq, and this contradiction proves JTSJ- Similarly, assuming (|T9"|) does not hold 
and defining 

zi = max {z e [-77, 0] | <j)\{z) = 0}, 

we conclude that z\ < and 4>'\( z i) < 0j so that 0a(-z) is decreasing in a 
neighborhood of zi, and this implies a contradiction to the definition of z\ 
and proves that (|79)l holds. This concludes the proof of the lemma. 



Since 



lim p c (X) 

A — >oo 



and p c (A) is a monotone function so that 

< p c (\) < e-i VA > 0, 
we conclude from lemma IT"!! that 

Lemma 15. In the excitable case f3 < 0, we have 

qe~2 

A ^W-« xJ(a;) * (A)<1 ' 

Let us note that since ^(A) < 1 implies that 1)52)1 doesn't hold, and since 
v — j, we can reformulate the previous lemma as a lower bound for the 
velocities of rotating waves in the excitable case. 

Lemma 16. In the excitable case (3 < 0, we have the following lower bound 
on the velocity of any rotating wave: 

2|/5|ei 1 

V > Tf \ ~- 



18 



GUY KATRIEL 



The following theorem shows that, in the excitable case and for sufficiently 
weak synaptic coupling, there are no rotating waves (so it implies part (II) (i) 
of theorem QJ. 

Theorem 17. In the excitable case /? < 0, if g £ (0,go), where 

n(c/3) 

90 — 77 \ ' 

max x( z R J{x) 

with fi(c, j3) defined by 162(1 . then there exists no rotating wave. 
PROOF: We shall show that if < g < go then 

(81) *(A) < 1 

for all A > 0, and thus that equation (|52jl cannot hold. By l|51|l. Ij81|l is 
equivalent to 

(82) a (tt) - a (-tt) < 2tt. 

We note that, by lemma El we already have l|81[) when (|77|l holds, hence we 
may assume 

(83) A<^/> c (A)maxJ(z). 

2\p\ xeR 

We define 

l" = & + 777Pc(A) max J (x) , 
and we note that (|83[) is equivalent to the statement that 

(84) fi > 0. 
Using 149(1 and lemma H~3l we have 

<p' x (z) = \[h{M*)) + gR\(z)w(M*))] 



< A 



1 - cos(4> x (x)) + (f3+ t?tPc(A) max J(xj) (1 + cos(0 A (z))) 
(85) = A[(/x + 1) + (m-1)cos(^(^))]. 

which implies (note that the integral below is well-defined because of (184(1 ) 

r m± < 2nX . 

J-n O* + 1) + (a*-1)«»(0a(«)) " 
Making the change of variables ip — 4>\{z) we obtain 

86 / 7 77 7-^— T7 7 — 7 < 2ttA. 

If we assume, by way of contradiction, that l|82l) does not hold, i.e. that 
0aM — <j>\{— 7r) > 27T, then, using (JSSJ, 



'^aC-t) (M + l) + cos(tp) 7 1) + (m- l)cos(yj) ^ 
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so together with ijHBt we obtain 

^ <2A, 

which is equivalent to 

1 1 - 4/3A 2 



> 



max ieS J(i) 2\p c {\) 

which contradicts g < go. This contradiction proves Ij82(l . concluding the 
proof of the theorem. 

We now proceed to prove that in the excitable case when the synaptic 
coupling is sufficiently large we have at least two rotating waves (see theorem 
1201 below). 

Lemma 18. In the excitable case (3 < 0, if there exists some Ao > with 

*(Aq) > 1, 

then there exist at least two solutions X\, A2 of l\5fy) with < A 2 < Ao < Ai, 
hence two rotating waves, with velocities satisfying 

1 1 1 
Vl = — < — < — = v 2 . 

Al Ao A2 

PROOF: By lemma IT21 we can choose A 2 < Ao so that ^(A^,) < 1. By lemma 
If 51 if we fix 

Ai = „, „, max J a: , 
1 2\f3\ xei v h 

then \&(A) < i- for all A > A' 1; and in particular it follows that A^ > Ao- We 
thus have 

X'2 <! Ao <C A^ 

with 

*(A 2 )<1, *(A )>1, *(A' 1 )<1. 

Thus by the intermediate value theorem, the equation l|52|) has a solution 
A 2 G (A 2 ,Ao) and a solution Ai e (Ao,Ai), corresponding to two rotating 
waves. 

The following lemma is valid for all values of (3: 
Lemma 19. Assume that A > satisfies the inequality 
(87) f C!0 (X)<gminJ(x), 

where f c p is defined by l\5y\) . Then 

#(A) > 1. 
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PROOF: By 1)5 our claim is equivalent to 

(88) ^(tt) - a (-tt) >2tt. 
We define 

?y = /3+ — p c (A) mm J(ar), 
and we note that (|87|l is equivalent to 

(89) V > ± 
Using 149|) and lemma IT31 we have 

i/ x (z) = \[h{M*)) + 9R&MM*))] 

> A[l - cos((/)x{x)) + (/3+J^p c (\)mmJ(x))(l + cos(<t>x(z))) 

(90) = A[(7 ? + l) + (r;-l)cos(0 A (z))], 
which implies 

M*)d* > z7rA . 



/-* (r? + l) + (r?-l)cos(^ A (2)) 
Making the change of variables <p — 4>\{z), we obtain 

(91) / 7 ^ , TT > 27rA - 

(r? + 1) + (?7 - 1) cos(^) " 

If we assume, by way of contradiction, that 1)88(1 does not hold, i.e. that 
</>a(7t) — </>a(— 7t) < 27r, then, using 

dy; . /" 27r dip 



/^(-tt) (»7 + !) + (V ~ 1) cos (v) Jo (?? + 1) + (r? - 1) cos(^) ^ 
so together with (|91l) we obtain 

— > 2A. 

This contradicts l|89[) . and this contradiction implies that (|88|) holds, complet- 
ing our proof. 

The following theorem implies part (II) (ii) of theorem ^ 
Theorem 20. In the excitable case (3 < 0, let 

(92) * = 

min ieR J{x) 

where f2(c, /?) is defined by \6ty) . Then when g > gi, there exist at least two 
rotating waves. In fact, we have a 'slow ' wave with velocity v s bounded from 
above by 

(93) v s < v cJgmin J(x)j 
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and a 'fast wave ' with velocity v / bounded from below by 

(94) v f >v c JgxamJ{x)), 

where v c p,v c ,(3 we the functions defined by |6'6)) . 

As a consequence of QU&p. \9^\j we have, for the slow wave 

p, „,<^£^I + o(^) 

mm xeR J{x) g Vg J / 
for the fast wave in the case c > 



as g oo, 



(96) v f> 2 \ 1 1 Vg + 0{1) asg^oo, 

y — 1 

and for the fast wave in the case c = 

(97) Vf > 2min J(x)g + o(^— ^ as <? — > oo. 

PROOF: g > gi and l|62[l imply the existence of A > satisfying l|87|) . hence by 
lemma [Tgl^fAl > 1, so that lemma IT8l implies the existence of two rotating 
waves. 

To prove l|93|) . i|94[l . we note that, assuming g > gi, the range of values of 
Ao for which (|87|1 holds is the interval 

Ac^femin J(x)) < A <X c ^(gxamJ(x)), 

where the functions A c «, A Cj( 3 are defined in section^ Thus, applying lemma 
dwith 

A = Ac.^gmin J(x)) - e, 

set 

where e > is arbitrarily small, we obtain the existence of a solution A £ of 
(|52[l with A £ > A c ,/3(<7min 2;e i{ J(x)) — e. Since e > is arbitrary, we have 
a solution A of (|52[1 with A > X c ^(gm.m x€ Tg_J(x)), hence a rotating wave 
with velocity v s satisfying i|93|l . Similarly applying lemma 1181 with Ao = 
A c g(g min^gR J( x )) + e > we obtain the existence of a wave with velocity Vf 
satisfying (|MJl . 

The estimates (ESJ-lEZIl follow from (^.(IMJ) and lemma ITUl 

We note that along with the upper bound l|95|l . we have a lower bound for 
the velocity of the slow wave, given by lemma [TBI 

Question 21. Derive an upper bound for the velocities of the fast waves (note 
that (Ml gives a lower bound). 

Question 22. Theoremsll7landl20lshow that several of the qualitative features 
that we saw explicitly in the case of uniform coupling (section^ remain valid 
in the general case. It is natural to ask whether more can be said, e.g., whether 



22 



GUY KATRIEL 



the following conjecture, or some weakened form of it, is true: for any J, there 
exists a value g C rit such that: 

(i) For g < g cr a there exist no travelling waves. 

(ii) For g = g cr n there exists a unique travelling wave. 

(iii) For g > g cr u there exist precisely two travelling waves. 

The next theorem deals with the oscillatory case j3 > 0, as well as the 
borderline case [3 = 0, and in particular proves part (I) of theorem ^ 

Theorem 23. // j3 > 7 there exists a rotating wave solution for any value of 
g > 0, with velocity v bounded from below by 

(98) v > v(q min J(x)), 

where v is the function defined by j71[ ), and the asymptotic formulas ifff61) , fl^ 
hold with Vf replaced by v. 

PROOF: If (3 > 0, then for any g > the equation 

/c,/3(A) = .9 mm J(x) 

has the unique solution A c ,p{g min^ga J{x)). Hence, any 

Ao > Xcpigvcan. J(x)) 

satisfies (|87f) . so that by lemma IT91 ^fAn) > 1. On the other hand for A 
sufficiently small we have, by lemma IT21 ^(A) < 1. Hence there exists a 
solution A € (0, Ao) of (|52p. Since Ao > A Ci ^(g min xe R J{x)) is arbitrary, we 
conclude that there exists a solution A < Xc.pig niin^gR J(x)) of l|52|l . Hence 
a rotating wave with velocity satisfying (|98ll . 

Question 24. Is it true that in the oscillatory case [3 > the rotating wave is 
always unique? We saw that this is the case when J = 1. 

Finally, we stress the important question of stability of the rotating waves, 
which remains open: 

Question 25. Investigate the question of stability of the rotating waves, i.e., 
do arbitrary solutions of (JSJ) , @ approach one of the rotating waves in large 
time? We conjecture that, at least under some restrictions on J, the rotating 
wave is stable in the case (3 > 0, while in the case [3 < the fast rotating 
wave is stable and the slow one is unstable. 
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